i/i 


AD-A145  752 
UNCLASSIFIED 


COMPARISON  OF  PREDICTED  AND  MEASURED  FRICTION  FOR  A 
TILT-PAD  THRUST  BERRI.  .  <U)  DAVID  W  TAVLOR  NAVAL  SHIP 
RESEARCH  AND  DEVELOPMENT  CENTER  BET.  T  L  DAUGHERTV 
JUL  84  DTNSRDC-84/825  F/G  20/11 


NL 


^  V  *. 


IM  12.5 


L25  1 1.4  HI  1.6 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BONCAU  or  STANOANOS-IBAS-A 


COMPARISON  OF  PREDICTED  AND  MEASURED  FRICTION  FOR  A  TILT-PAD  THRUST 
BEARING  WITH  CROWN  AND  ROUGHNESS  UNDER  HYDRODYNAMIC  AND 

MIXED  LUBRICATION  DTNSRDC-84/025 


AD- A 145  752 


DAVID  W.  TAYLOR  NAVAL  SHIP 
RESEARCH  AND  DEVELOPMENT  CENTER 

Bathaada,  Maryland  20094 


COMPARISON  OF  PREDICTED  AND  MEASURED  FRICTION  FOR 
A  TILT-PAD  THRUST  BEARING  WITH  CROWN  AND 
ROUGHNESS  UNDER  HYDRODYNAMIC  AND 
MIXED  LUBRICATION 


by 

Thomas  L.  Daugherty 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED 

) 


SHIP  MATERIALS  ENGINEERING  DEPARTMENT 
RESEARCH  AND  DEVELOPMENT  REPORT 


DT1C 


( 

jL 


MAJOR  DTNSRDC  ORGANIZATIONAL  COMPONENTS 


NDW-DTNSRDC  3960/43b  («•».  2-801 


SECURITY  CLASSIFICATION  of  This  PAGE  ('IFh#n  Dote  Entered) 


REPORT  DOCUMENTATION  PAGE 


1  REPORT  NUMBER 


84 / 022 


4  TITLE  (and  Subtitle)  S.  TYPE 

COMPARISON  OF  PREDICTED  AND  MEASURED  FRICTION  FOR  Final 
A  TILT-PAD  THRUST  BEARING  WITH  CROWN  AND 
ROUGHNESS  UNDER  HYDRODYNAMIC  AND 
MIXED  LUBRICATION 


7.  authorc*; 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


3  RECIPIENT’S  CATALOG  NUMBER 


S.  TYPE  OF  REPORT  ft  PERIOD  COVERED 


6.  PERFORMING  ORG.  REPORT  NUMBER 


Thomas  L.  Daugherty 


a.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 


S  CONTRACT  OR  GRANT  NUMBER/*,) 


10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  ft  WORK  UNIT  NUMBERS 


It.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 


David  W.  Taylor  Naval  Ship  Research  &  Development  Program  Element  61152N 
Center,  Bethesda,  MD  20084  Task  Area  ZR0240301 

Work  Unit  2832-121 


12.  REPORT  DATE 

July  1984 


13.  NUMBER  OF  PAGES 

52 


14.  MONITORING  AGENCY  NAME  ft  ADDRESS/!/  dllttnnl  from  Controlling  Olllet)  IS.  SECURITY  CLASS,  (ot  thla  rmport) 

UNCLASSIFIED 


ISa.  DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 


16.  DISTRIBUTION  STATEMENT  (ol  thlt  Rmport) 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED. 


17.  DISTRIBUTION  STATEMENT  (of  the  mbetrect  entered  in  Stock  20,  It  different  from  Report) 


19-  KEY  WORDS  (Continue  on  reveree  tide  If  neceeeery  end  Identity  by  Woe*  number) 


Bearings 
Hydrodynamic : 
Lubrication 


Mixed  Lubrication 
Friction 
Thrust  Bearings 


20.  ABSTRACT  (Continue  on  reveree  elde  It  neceeeery  end  Identity  by  block  number) 

4  In  a  controlled  experiment  on  a  tilt-pad  thrust  bearing 
having  measured  surface  crown  and  longitudinal  roughness, 
the  frictional  results  were  compared  with  results  predic¬ 
ted  by  a  physico-mathematical  model  of  the  bearing.  The 
bearing  operated  hydrodynaraically  using  2190  TEP  oil  at  ^ 

(Continued  on  reverse  side) 


DD 


FORM 
JAN  73 


EDITION  OF  I  NOV  «S  IS  OBSOLETE 
S/N  0  102-0  14-  660  1 


(UNCLASSIFIED) 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whan  Dad  > irimd) 


_ UNCLASSIFIED _ 

__'_UW|TY  CLASSIFICATION  OF  this  PAGE'lthen  Data  Entond) 


Block  20  continued  ^ 

•^“ambient  and  150*F  bulk  temperature  under  loads  approaching 
500  psiV^The  curve  of  the  coefficient  of  friction  as  a 
function  of  the  parameter  pU/PB  derived  from  the  experiment 
was  parallel  to  the  curve  predicted  by  the  model  under 
hydrodynamic  conditions  but  the  values  were  higher  than, 
those  predicted.  This  difference  was  attributed  primarily 
to  end  leakage,  which  was  not  considered  in  the  model. 

A  low  viscosity  lubricant,  MIL-H-5606 ,  was  used  in  experi¬ 
ments  in  the  transition  between  hydrodynamic  and  mixed 
lubrication  regimes.  However,  this  fluid  exhibited  non- 
Newtonian  behavior  thereby  complicating  comparisons  with 
model  predictions. 

Because  so  little  experimental  data  are  available  on 
the  tribological  interaction  occurring  with  non-Newtonian 
fluids,  additional  experiments  should  be  conducted  so  that 
a  more  sophisticated  model  can  be  developed  of  slider 
bearings  such  as  the  tilt-pad  thrust  bearing  operating  in 
the  mixed  lubrication  regime. 

Some  areas  in  need  of  additional  work  were  outlined, 
and  experiments  are  needed  to  confirm  some  of  the  assumptions 
used  in  model  analyses. 


_ UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  THIS  PAGEfRTian  Dafa  Enttrud) 


TABLE  OF  CONTENTS 


LIST  OF  FIGURES . iii 


LIST  OF  TABLES 


NOMENCLATURE  .  vi 


ABSTRACT 


ADMINISTRATIVE  INFORMATION 


BACKGROUND 


OBJECTIVE 


APPROACH 


THEORETICAL  BASIS  FOR  THE  MODEL 


EXPERIMENTAL  MATERIALS  .  11 

TEST  MACHINE . 11 

DESCRIPTION  OF  THE  TEST  BEARINGS . 15 

TOPOGRAPHICAL  ANALYSIS  OF  THE  TEST  BEARINGS  .  18 

COMPUTER  ANALYSIS  .  26 

EXPERIMENTAL  RESULTS  .  33 


DISCUSSION 


CONCLUSIONS . 40 

RECOMMENDATIONS  .  40 

ACKNOWLEDGEMENTS . 41 

REFERENCES . 43 


LIST  OF  FIGURES 


1  -  Stribeck  Curve  and  Lubrication  Regimes 


2  -  Curved  Slider  Bearing  with  Surface 
Roughness  . 


3  -  Crowned  Bearing  Surface  with  Regions 

of  Hydrodynamic  and  Mixed  Lubrication  ...  .  9 


s 


Page 


4  -  Test  Machine  Used  to  Conduct  Friction 

Measurements . 12 

5  -  Test  Machine  Schematic .  13 

6  -  Thrust  Runner  Showing  Octagonal  Bearing 

and  Mating  Surface  for  Top  Set  of  Bearings .  14 

7  -  Three-Pad  Thrust  Bearing  Mounted  on  Leveling 

Washer .  15 

8  -  Front  and  Back  Surfaces  of  Tilt-Pad  Bearing .  16 

9  -  Thrust  Bearing  Assembly  Drawing  .  17 

10  -  Grids  for  Topographical  Measurements  of 

Thrust  Pads  and  Runner  Surfaces .  19 

11  -  Thrust  Pad  Talysurf  Traces  at  10,000X 

Vertical  and  2X  Horizontal  Magnification 

for  Bearing  T-2,  Outside  Axial  Trace  .  20 

12  -  Thrust  Runner  Talysurf  Traces  at  10,000X 

Vertical  and  2X  Horizontal  Magnification  for 

Thrust  Runner  A  at  270°  Position .  20 

13  -  Characterization  of  Bearing  Surface  Traces 

on  Talysurf  10 . 21 

14  -  Typical  Computer  Output  Under  Hydrodynamic 

Lubrication  Conditions . 27 

15  -  Typical  Computer  Output  Under  Mixed 

Lubrication  Conditions  .  . . 28 

16  -  Pressure  Distribution  Under  Hydrodynamic 

Lubrication  Conditions  Allowing  Negative 

Pressures .  29 

17  -  Pressure  Distribution  Under  Hydrodynamic 

Lubrication  Conditions  with  Reynolds-Swif t- 

Stieber  Boundary  Conditions  .  29 


18  -  Pressure  Distribution  Under  Mixed  Lubrication 

Conditions  with  Reynolds-Swif t-Stieber  Boundary 
Conditions . . . 

19  -  Pressure  Distribution  Under  Mixed  Lubrication 

Conditions  Allowing  Negative  Pressures  .... 


iv 


Theoretical  Friction  Coefficients  Versus 
Parameter  yU/PB  Allowing  Negative  Pressure 
and  Applying  Reynolds-Swif t-Stieber 
Boundary  Conditions  . 

Theoretical  Minimum  Film  Thickness  Versus 
Parameter  yU/PB  Allowing  Negative  Pressure 
and  Applying  Reynolds-Swif t-Stieber 
Boundary  Conditions  . 

Viscosity-Temperature  Relationship  for 
MS  2190  TEP  Oil  and  MIL-H-5606  Hydraulic 
Fluid  . 

Comparison  of  Measured  and  Computed  Friction 
Coefficient  for  Crowned  Tilt-Pad  Thrust 
Bearing  Having  Longitudinal  Roughness  .  .  .  . 

LIST  OF  TABLES 

Roughness  and  Crown  Measurements  for  Bottom 
Set  of  Thrust  Pads  on  a  Tilt-Pad  Bearing  .  . 

Roughness  and  Crown  Measurements  for  Top 
Set  of  Thrust  Pads  on  a  Tilt-Pad  Bearing  .  . 

Roughness  and  Crown  of  Thrust  Collars  .  .  .  . 

Maximum  Recorded  2c  Values  for  the  Thrust 
Runner  and  Bearing  Surfaces  . 

Computer  Input  Based  on  Topography  of 
Surfaces  .  ....... 

Data  from  Run  1  -  2190  TEP  at  Ambient 
Temperature  . 

Data  from  Run  2  -  2190  TEP  at  150°F  . 

Data  from  Run  3  -  MIL-H-5606  Fluid  at 


NOMENCLATURE 


Ratio  of  inlet  to  outlet  film  thickness  (hj/f^) 

Bearing  width  in  the  direction  of  motion 

Half  total  range  of  random  film  thickness  variable 

Constant  of  integration 

Center  of  pressure 

First  derivative  with  respect  to  x 


E  (  )  Expected  Value  of  (  ) 


Coefficient  of  friction 


Total  applied  thrust  load 
Normal  applied  force 


Runner  friction  force 


Force  to  shear 


Dummy  function  of  x 


Smooth  fluid  film  thickness  as  function  of  x 


Smooth  film  thickness 

Inlet  film  thickness  to  bearing 

Outlet  film  thickness  to  bearing 

Limiting  value  of  h  used  to  prevent  In  (h-h*)  from  being  unbounded 
Minimum  film  thickness  for  fixed  inclined  slider 
Randomly  distributed  roughness  height 
Film  thickness  including  roughness  deviations 


Bearing  pad  crown  height 


A  constant 


33 


Friction  coefficient  factor 


Pressure  coefficient  factor 


ABSTRACT 


t 


In  a  controlled  experiment  on  a  tilt-pad  thrust  bearing 
having  measured  surface  crown  and  longitudinal  roughness, 
the  frictional  results  were  compared  with  results  predic¬ 
ted  by  a  physico-mathematical  model  of  the  bearing.  The 
bearing  operated  hydrodynamically  using  2190  TEP  oil  at 
ambient  and  150°F  bulk  temperature  under  loads  approaching 
500  psi.  The  curve  of  the  coefficient  of  friction  as  a 
function  of  the  parameter  PU/PB  derived  from  the  experiment 
was  parallel  to  the  curve  predicted  by  the  model  under 
hydrodynamic  conditions  but  the  values  were  higher  than 
those  predicted.  This  difference  was  attributed  primarily 
to  end  leakage,  which  was  not  considered  in  the  model. 

A  low  viscosity  lubricant,  MIL-H-5606 ,  was  used  in  experi¬ 
ments  in  the  transition  between  hydrodynamic  and  mixed 
lubrication  regimes.  However,  this  fluid  exhibited  non- 
Newtonian  behavior  thereby  complicating  comparisons  with 
model  predictions. 

Because  so  little  experimental  data  are  available  on 
the  tribological  interaction  occurring  with  non-Newtonian 
fluids,  additional  experiments  should  be  conducted  so  that 
a  more  sophisticated  model  can  be  developed  of  slider 
bearings  such  as  the  tilt-pad  thrust  bearing  operating  in 
the  mixed  lubrication  regime. 

Some  areas  in  need  of  additional  work  were  outlined, 
and  experiments  are  needed  to  confirm  some  of  the  assumptions 
used  in  model  analyses. 
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ADMINISTRATIVE  INFORMATION 

This  report  covers  work  conducted  for  and  funded  by  David  Taylor  Naval  Ship 
Research  and  Development  Center's  Independent  Research  project  entitled  "Mixed 
Lubrication  Model."  The  study  was  carried  out  under  Program  Element  61152N,  Task 
Area  ZR0240301,  and  Work  Unit  2832-121. 

BACKGROUND 

A  wide  variety  of  sliding  surface  components  used  in  shipboard  machinery  are 
designed  to  operate  with  a  film  of  lubricant  between  the  mating  surfaces.  The  bear¬ 
ings  on  the  main  shafting  are  examples.  Inboard  bearings  usually  consist  of  an  oil 
lubricated  main  thrust  bearing  and  several  oil  lubricated  lineshaft  journal  bearings, 
whereas  outboard  bearings  in  the  U.S.  Navy  are  lubricated  with  seawater  and  are 
stave  type  or  full  molded  sterntube  and  strut  bearings.  All  these  bearings  are  of 
the  hydrodynamic  type.  That  is,  lubricant  film  is  established  by  the  relative 
sliding  velocity  between  the  bearing  surface  and  the  mating  surface.  The  lubricant 
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film  that  is  developed  forces  the  two  surfaces  apart  and  thereby  eliminates  both 
surface  contact  and  wear.  However,  these  bearings  must  operate  part  of  their  life 
under  conditions  where  there  is  no  full  lubricant  film.  Starts,  stops,  transient 
operations,  and  operation  on  jacking  gear  are  examples  of  some  of  those  conditions. 
If  lubricant  films  are  not  thick  enough  to  keep  the  surfaces  apart,  surface 
irregularities  interact  and  wear  results. 

The  frictional  performance  of  such  bearings  can  be  described  by  use  of  the 
Stribeck  curve-*-*  shown  in  Figure  1.  Three  lubrication  regimes  are  shown —  bounda",r 
mixed,  and  hydrodynamic.  The  curve  shows  that  friction  is  dependent  upon  lubricant 
viscosity  (p),  the  relative  velocity  of  the  mating  surfaces  (U),  and  the  normal 
force  (Ffj). 


BOUNDARY  MIXED  HYDRODYNAMIC 


VISCOSITY,  p  x  VELOCITY.  U 
LOAD,  Fn 

Figure  1  -  Stribeck  Curve  and  Lubrication  Regimes 


lf|  —  ■■  '  . ■■■■■■  . 

References  are  listed  on  page  43. 
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In  the  hydrodynamic  region,  the  rigid  surfaces  are  separated  by  a  continuous 
lubricant  film,  whose  thickness,  h,  is  much  larger  than  the  combined  surface 
roughness,  Rq.  Friction  is  primarily  due  to  shearing  within  the  lubricant.  This 
is  the  desirable  region  of  operation  for  most  bearings  because  there  is  no  contact 
and  thus  no  wear. 

There  is  contact  between  asperities  if  the  viscosity  of  the  lubricant  or  the 
relative  sliding  velocity  is  reduced  or  if  the  load  is  increased  so  that  the  film 
thickness  approaches  the  combined  surface  roughness,  Rq.  Then  the  load  is  carried 
by  lubricant  film  pressure  produced  by  hydrodynamic  effects  as  well  as  by  asperity 
contact.  Under  these  mixed  lubrication  conditions,  the  magnitude  of  load  carried 
by  asperities  increases  with  lower  values  of  pU/Fjg.  Friction  and  wear  also  increase 

If  pU/Fn  is  such  that  the  film  thickness  approaches  zero,  almost  all  the  load 
is  carried  by  the  asperity  contacts.  Hydrodynamic  influence  is  negligible  and  wear 
is  governed  more  by  the  physico-chemical  properties  of  the  lubricant  and  the  mating 
surface  and  their  interaction. 

Design  charts  have  been  developed  for  frequently  used  bearings  such  as  the 
tilt-pad  slider  bearing. 2  Operating  characteristics  of  the  tilt-pad  bearing  de¬ 
pend  on  the  pivot  position  and  the  profile  of  the  pad  surface  and  can  be  optimized 
by  the  proper  selection  of  pivot  and  profile.  The  tilt-pad  thrust  bearing  the 
U.S.  Navy  uses  for  the  main  propulsion  shafting  is  pivoted  in  the  center  of  the 
pad, thereby  allowing  good  (not  optimum)  operation  in  either  direction.  According 
to  simple  theory  using  the  Reynolds  equation,  a  centrally  pivoted  pad  with  a  flat 
surface  has  no  load  capacity. 2  However,  experience  has  shown  that  centrally 
pivoted  pads  perform  reasonably  successfully.  Load  capacity  of  a  centrally 
pivoted  tilt-pad  bearing  can  be  predicted  adequately  by  considering  effects  of 
pad  crown, 3  of  variable  lubricant  viscosity  as  the  lubricant  passes  between  the 
bearing  and  the  mating  surface,  and  of  variable  lubricant  density. 

Crowning  of  the  pad  surface  has  a  strong  influence  on  improving  the  performance 
of  the  centrally  pivoted  tilt-pad  bearing.  Some  degree  of  crowning  is  expected  on 
new  bearings,  but  most  results  from  deformations  caused  by  the  oil  film  pressure 
distribution  and  thermal  gradient  across  the  pad  thickness  during  operation. 

However,  an  excessive  crown  may  produce  poor  performance. 
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Although  the  frictional  performance  of  tilt-pad  bearings  operated  in  the  hydro- 
dynamic  lubrication  regime  can  be  predicted  by  theoretical  consideration,  the 
literature  contains  little  experimental  evidence  of  performance  in  the  transition 
region  between  hydrodynamic  and  mixed  lubrication  conditions.  Thus,  the  goal  of 
the  project  was  to  develop  a  model  for  sliding  surface  bearings  that  would  extend 
the  hydrodynamic  theory  continuously  into  the  mixed  lubrication  region.  This 
report  compares  the  frictional  results  from  a  controlled  experiment  on  a  tilt-pad 
thrust  bearing  having  measured  surface  and  longitudinal  roughness  with  the 
predictive  results  from  a  model  of  the  bearing. 

APPROACH 

A  set  of  centrally  pivoted  tilt-pad  thrust  bearings  was  obtained  commercially. 
The  pads  were  crowned  to  simulate  thermal  and  elastic  pad  distortions  and  roughened 
in  a  longitudinal  direction  to  simulate  expected  roughness  patterns  generated  during 
operation.  The  topography  of  each  pad  and  mating  runner  surface  was  measured  and 
the  surface  characterized  by  crown  height  and  a  roughness  magnitude.  These  values 
and  bearing  and  material  description,  speed,  and  lubricant  characteristics  were  used 
in  a  computer  program  of  a  model  of  the  bearing.  The  theory  and  the  model  are 
described  below. 

The  predicted  frictional  behavior  of  each  pad  and  of  the  complete  bearing  were 
obtained  from  the  computer  program.  The  resulting  friction  coefficient  was  plotted 
against  the  parameter  in  pU/PB  to  form  a  plot  similar  to  the  Stribeck  curve  shown  in 
Figure  1.  A  friction  curve  resulted  showing  a  prediction  of  performance  in  the 
hydrodynamic  region  of  operation  and  into  the  mixed  lubrication  region.  Other 
meaningful  outputs,  such  as  the  minimum  film  thickness,  inlet  to  exit  film  thickness 
ratio  and  location  of  mixed  lubrication  zone,  were  also  generated  by  the  computer 
program. 

Next  a  series  of  experiments  was  conducted  on  the  set  of  modified  thrust 
bearings  over  a  range  of  pU/PB  values  and  the  frictional  torque  recorded.  The 
coefficients  of  friction  were  computed  for  the  values  of  pU/PB  and  compared  with 
the  predicted  results. 


THEORETICAL  BASIS  FOR  THE  MODEL 


The  model  for  mixed  and  hydrodynamic  lubrication  of  a  crowned  tilt-pad  thrust 
bearing  with  longitudinal  roughness,  which  was  developed  under  this  project,  is  an 
extension  of  an  earlier  analysis  of  a  fixed  inclined  bearing  with  longitudinal 
roughness.^  The  bearing  pad  is  free  to  pivot  until  a  moment  equilibrium  position 
is  attained.  The  model  geometry  is  shown  in  Figure  2.  The  pad  surface  is  crowned 
with  the  maximum  pad  height  at  the  center.  The  film  thickness  due  to  the  crowning 
the  angle  of  inclination,  the  ratio  of  inlet  to  outlet  film  thickness,  and  the  pad 
geometry  can  be  expressed  as: 

h(x)=h2  {a  +  (1-a)  [(1  +  4  Hc/mB)(x/B)  -  (4  Hc/mB)g)2]} 

The  expression  is  for  a  smooth  surface  and  is  a  function  of  the  x  coordinate. 
Another  term  must  be  added  when  roughness  is  being  considered.  A  film  height 
would  then  consist  of: 

H  =  h(x)  +  hs 

where  hs  is  the  randomly  distributed  roughness  term. 


(NOTE:  IF  BOTH  SURFACES  ARE  ROUGH,  2C  IS 
THE  COMBINED  ROUGHNESS  HEIGHT) 

Figure  2  -  Curved  Slider  Bearing  with 
Surface  Roughness  (from  Daugherty  and  Pan^) 


Experimental  evidence  indicates  that  many  engineering  surfaces  have  a  distribution 
of  random  heights  that  is  very  nearly  Gaussian,  at  least  within  three  standard 
deviations.  We  used  Christensen's  polynomial  distribution  function  to  approximate 
the  Gaussian  distribution. ^ 

f  (hg)  =  15  (c2  -  hg2)3,  -  c<hs<c 

32c7 

=  0  elsewhere 

where  2c  is  the  domain  of  the  function.  The  function  becomes  zero  at  c=3o, 
where  o  is  the  standard  deviation.  The  above  equation  will  be  considered  a  descrip¬ 
tion  of  the  combined  roughness  when  both  the  bearing  and  runner  surfaces  are  rough 
The  following  assumptions  were  made: 

•  Lubricant  viscosity  is  constant. 

•  Lubricant  density  is  constant. 

•  Reynolds'  equation  holds  for  hydrodynamic  and  mixed  lubrication 
conditions. 

•  There  is  no  side  leakage. 

•  No  effects  are  caused  by  material/lubricant  interactions  of  a  chemical- 
physical  nature. 

•  Centrifugal  effects  are  negligible. 

•  Lubricant  is  Newtonian. 

•  Film  pressure  does  not  vary  across  film. 

The  basic  form  of  Reynolds'  equation  under  steady  state  conditions  is 

d_  E  (H3)  d£  =  6UU  d  E(H) 

dx  dx  dx 

Where  E(  )  is  the  expected  value  of  (  ) ,  and  H  =  h(x)  +  hs.  The  quantity  p  is  the 
mean  hydrodynamic  pressure  at  any  location  x,  U  is  the  slider  speed,  and  p  is  the 
lubricant  viscosity.  The  expectancy  operator  is  defined  as: 

E[g(x) ]  =  J  g(x)f(hs)dhs, 

where  the  upper  limit  of  integration  is  always  c,  and  the  lower  limit  -c  or  -h,  de¬ 
pending  upon  whether  the  film  thickness  is  greater  or  less  than  c.  This  corresponds 
to  regions  of  pure  hydrodynamic  and  mixed  lubrication.  In  regimes  of  pure  hydro- 
dynamic  lubrication,  it  can  be  shown  that: 


E(K)  =  K,  a  constant 


E(H)  =  h, 

and  E(H3)  =  h3  +  h  c2/3 

Under  mixed  lubrication  conditions  these  parameters  are  much  more  complicated: 


E(K)  =  K/32c7)  16c7  +  35c8h  -  35c4h3  +  21c2h3  -  5h7 


E(H)  =  35/ 32c7)  £8  +  16c7h  +  J_c8h2  -  lc*h4  +  _lC2h8  -  _JT*8 

8  35  2  4  10  56 


E(H3)=35/32c7  lc10+  16c9h+3c8h2+16c7h3+lcM-  lc4h6+  3  c2h6-  1  hlO 

40  105  8  35  4  20  280  840 


The  pressure  distribution  within  the  lubricant  film  can  be  determined  by 
integrating  Reynolds'  equation.  The  first  integration  yields: 


d£  =  6uU  |(E(H)/E(H3)  +  C/E(H3)J 


where  C  is  a  constant  of  integration. 
The  second  integration  yields: 


P  =  6vU  /  (E(H)/E(H3))dx 
0 


x  /  \ 
+  Cf  1  \ 

0  \E(H3)/ 


The  boundary  conditions  are  p(o)  =  p(B)  =  0.  The  constant  of  integration  is: 


B  B 

C  -  -  /  (E(H)/E(H3))dx/  / 
0  0 


/// 1  \ 

o  y e(h3) J 


In  the  mixed  lubrication  region,  the  metallic  contact  load  is  assumed  equal  to 
the  product  of  the  yield  stress  of  the  weaker  bearing  material  and  the  fraction  of 
area,  per  unit  bearing  width,  over  which  such  contact  takes  place.  The  load 
capacity  per  unit  width  can  therefore  be  expressed  as: 


1]  lo  B 

w  -  /  pdx  +J  pdx  +  J  pdx 
0  li  12 


I 


where  p  is  the  ordinary  hydrodynamic  pressure  upstream  and  downstream  of  the  area 
of  asperity  interactions  and  p  is  the  pressure  in  the  mixed  lubrication  zone.  These 
regions  are  shown  visually  in  Figure  3.  In  the  first  region  from  x=0  to  x=lj,  the 
film  thickness  exceeds  the  critical  film  thickness  and  converges  in  the  direction  of 
motion.  Significant  load  carrying  hydrodynamic  pressure  is  generated  in  the  lubri¬ 
cant  film.  In  the  region  from  x=l^  to  x=l2  the  film  thickness  is  less  than  the 
critical  film  thickness,  producing  mixed  lubrication.  Load  support  is  generated  by 
a  combination  of  asperity  contact  and  hydrodynamic  pressure  between  asperity 


interactions.  In  this  region  the  smooth  film  thickness  goes  through  a  converging 


and  then  a  diverging  region.  In  the  remaining  region  from  x=l2  to  x=B  the  film 
thickness  again  exceeds  the  critical  film  thickness,  thereby  signifying  that  there 
are  no  asperity  interactions.  The  pressure  produced  in  this  region  results  from 
hydrodynamic  action  only.  The  film  thickness  diverges  in  the  direction  of  motion; 
therefore  no  load  support  is  generated  in  this  region. 

Load  carried  by  the  second  region  is  supported  in  part  by  the  hydrodynamic 
pressures  generated  between  the  asperities  that  are  in  contact  and  in  part  by  the 
asperity  loads  themselves  as  follows: 


\2  \2 

f  pdx  =  J  r(x)pdx  +  py  J  (l-r)dx 

U  Ii  Ii 


where  Py  is  the  yield  stress  of  the  weaker  bearing  material  and  r  is  the  fraction  of 
bearing  area  formed  by  the  valleys  between  asperity  elements 


1  16c7  +  35c&h  -  35c*h3  +  21c2h3  -  5h7 j 

32c7 


per  unit  bearing  width. 

The  center  of  pressure  can  be  obtained  by  the  following  formula: 


B  B 

f pxdx/  f  pdx 

0  0 


The  center  of  pressure  is  used  to  obtain  moment  balance  of  the  thrust  pad  about  its 
center  of  pressure  by  'equating  the  product  of  the  load  and  the  center  of  pressure 
distance  from  the  pivot  point  to  the  product  of  the  friction  on  the  pad  surface  and 
its  distance  to  the  center  of  pivot. 

The  shear  stress  acting  on  the  slider  surface  from  hydrodynamic  action  is  given 

by: 


TOSS* 


.'VW 


n.n.iw'.rvxn 
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x(x)  *  pUE(l/H)  +  2  l^dx/  E^H' 

The  total  friction  force  can  be  determined  in  a  similar  manner  to  the  load  capacity. 


That  is. 


ll  I2  ^  B 

>  f  Tdx  +J  xdx  +/  Tdx 
0  ll  12 
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where  f  idx  =j  rxdx  +  s  f  (1-r) 
ll  ll  ll 


where  s  is  the  shear  stress  of  the  weaker  bearing  material. 

The  expectation  value  E(l/H)  must  be  evaluated  for  the  hydrodynamic  and 
mixed  regions  to  determine  the  friction  force  - 


j/A  =  15  (c2  -  h2)3  In  (h  +  c)  +  _2  ch(15h4  -  40c2h2 

l H j  32c7  (h  -  c)  15 


+  33c4) 


for  the  hydrodynamic  region  and: 


E/l\  =  35  (c2  -  h2)3  In  (h  +  c  )|+  J3h3(h  +  c)3  -  l^(h  +  c)6 

W  32c7  (h  -  h*)  3  6 


-  4(c2  -  h2)h(h  +  c)3  +  6(c2  -  h2)2h(h  +  c)  +  6(c2  -  h2)h2(c  +  h)2  -  3h2(c  +  h)4 

+  6h(h  +  c)3  +  2  (c^  “  h2)(h  +  c)4  -  3(c^  ”  h2)(h  +  c)2] 

5  4  2 

for  the  mixed  region. 
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Figure  3  -  Crowned  Bearing  Surface  with  Regions 
of  Hydrodynamic  and  Mixed  Lubrication 


The  term  h*  is  needed  to  prevent  the  logarithmic  term  from  being  unbounded;  h*  was 
set  eaual  to  h-mB/1000.  Numerical  analysis  showed  that  the  resulting  solutions  were 
not  particularly  sensitive  to  this  assumption. 

The  coefficient  of  friction  may  be  determined  by  dividing  the  friction  force 
by  the  applied  load.  The  friction  force  is  approximated  as  the  sum  of  the  friction 
forces  (as  determined  in  the  foregoing  paragraphs)  and  the  component  of  load  acting 
in  the  plane  of  rotation  or: 

Fr  =  Fg+mW  where  Fr  is  the  force  on  the  runner. 

The  friction  coefficient,  f,  is  then: 

f  =  Ir  m  Is  +  m 

U  W 

The  inputs  to  the  computer  program  include  the  rotational  speed,  the  pad  length, 
the  viscosity  of  the  lubricant,  the  mean  bearing  surface  diameter,  the  roughness 
asperity  height,  the  bearing  material  yield  strength,  the  ratio  of  the  shear  stress 
to  yield  stress,  the  crown  height,  and  the  pivot  coordinates. 

The  computer  program  is  written  in  BASIC  for  an  interactive  system.  It  calcu¬ 
lates  a  reasonable  first  approximation  of  h]/h2.  Reynolds'  equation  is  then 
successively  integrated  to  obtain  the  friction  and  pressure  moments.  If  the  moment 
is  greater  than  the  product  of  the  load  and  one-tenth  of  one  percent  of  the  bearing 
length,  then  a  new  value  of  h^/t^,  calculated  from  a  “hunting"  routine,  is  inserted, 
and  the  process  is  repeated  until  the  moment  balance  criterion  is  met. 

When  the  moment  balance  criterion  has  been  met,  the  operator  receives  a  list  of 
pertinent  bearing  parameters  such  as  load  and  friction  forces  broken  down  into 
contributions  from  pure  hydrodynamic,  mixed,  and  asperity  contact  forces;  shear  and 
pressure  moments;  the  center  of  pressure;  coefficient  of  friction;  value  of  uU/PB; 
hi/h2  ratio;  and  the  minimum  film  thickness  and  its  location.  The  operator  can 
choose  additional  conditions  and  rerun. 

The  initial  computer  program  was  developed  by  Dr.  B.  Carson  under  Navy  Contract 
N00600-78-D-1302  (HW-22)  with  CADCOM,  a  division  of  Mantech  International 
Corporation.  That  program  was  modified  to  improve  the  iteration  scheme  used  (which 
reduced  computer  time  and  improved  the  search  portion);  to  obtain  a  pressure  plot 
of  the  film  within  the  bearing;  and  to  determine  effects  of  cavitation  in  the 
diverging  region  of  the  film.  In  the  original  program,  significant  negative 


pressures  were  predicted  using  Reynolds'  equation  and  boundary  conditions  of 
p(0)  =  p(B)  =0.  A  diverging  film  was  formed  when  the  trailing  edge  was  not  the 
location  of  minimum  film  thickness.  It  is  generally  believed  that  a  subambient 
pressure  condition  in  the  fluid  cannot  be  sustained ; 6 instead,  the  film  ruptures 
or  cavitates  somewhere  so  that  the  ambient  condition  will  be  reached  internally 
with  a  vanishing  pressure  gradient.  This  is  the  Reynolds-Swif t-Stieber  condition 
of  film  rupture  or  cavitation  and  is  stated  as: 

p(xc)  =  0 

dp  =  0 

dx  x 

c 

0  <  xc  <  B 

where  xc  is  the  location  of  the  cavitation  or  rupture  boundary. 

EXPERIMENTAL  MATERIALS 

TEST  MACHINE 

Figure  4a  is  a  photograph  of  the  machine  we  used  to  measure  friction  on  crowned 
tilt-pad  bearings  having  a  longitudinal  roughness  pattern.  The  machine  has  six 
main  parts — hydraulic  motor,  reduction  gear,  torque  sensor,  oil  tank,  load  blocks, 
and  heater.  The  test  bearing  (Figure  4b)  is  located  in  the  oil  tank  section.  A 
schematic  of  the  test  machine  is  shown  in  Figure  5.  The  shaft  is  vertical  and  is 
positioned  with  ball  bearings.  An  octagonal  bearing  with  a  series  of  linear  bear¬ 
ings  on  each  face  allows  the  runner  to  move  axially  while  transmitting  torque 
(Figure  6).  Two  sets  of  three  tilt-pad  bearings  are  mounted  on  a  metal  ring  and  are 
submerged  in  the  lubricant.  The  oil  reservoir  contains  7-1/2  quarts  of  lubricant 
and  the  frame  is  made  of  plexiglass  to  allow  viewing  during  operation.  Two  steel 
plates  are  attached  to  the  runner  as  mating  surfaces  to  the  bearings.  The  lubricant 
can  be  heated  by  circulating  hot  water  through  coils. 


1  HYDRAULIC  MOTOR 

2  REDUCTION  GEAR 

3  TORQUE  SENSOR 


4  OIL  TANK 

6  LOAD  PISTON  BLOCKS 

5  HEATER 


FIGURE  4a  -  TEST  MACHINE 


1  HEATING  COIL  3  TEST  BEARING 

2  RUNNER  PLATES  4  THERMOCOUPLE 


FIGURE  4b  -  TEST  BEARING  IN  OIL  TANK 


Figure  4  -  Test  Machine  Used  To  Conduct 
Friction  Measurements 
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Figure  5  -  Test  Machine  Schematic 
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Figure  6  -  Thrust  Runner  Showing  Octagonal  Bearing 
and  Mating  Surface  for  Top  Set  of  Bearings 


The  rotational  speed  is  set  at  a  nominal  35  rpm.  A  hydraulic  load  is  applied 
to  the  upper  set  of  bearings  and  transmitted  through  film  to  the  runner,  back 
through  the  film  of  the  lower  bearing  set  and  through  the  lower  set  to  the  test 
machine  frame.  The  thrust  load  is  carried  by  only  the  two  sets  of  test  bearings. 
This  allows  accurate  measurements  of  friction  in  test  bearings  because  support 
bearings  cause  little  torque.  The  parameters  measured  are  speed,  torque,  load,  and 
bulk  oil  temperatures.  The  thrust  runners  are  of  1040  steel  (50  Rc  Hardness).  Thev 
are  machined  to  produce  circumferential  machining  marks,  which  represent  the  desired 
longitudinal  roughness  pattern.  Roughness  measurements  are  described  in  the  section 
on  topographic  analysis  below.  The  thrust  runner  is  shown  in  Figure  h. 

Load  is  applied  by  supplying  hydraulic  pressure  to  three  pistons  located  in  the 
top  load  block.  Each  piston  is  located  directly  over  the  thrust  bearings.  A  load 
cell  in  each  piston  monitors  the  applied  load. 
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DESCRIPTION  OF  THE  TEST  BEARINGS 


The  bearing  we  used  in  these  experiments  is  a  commercially  available,  three 
shoe  element  tilt-pad  thrust  bearing  made  by  Kingsbury  Machine  Works,  Inc.,  as 
bearing  number  NN-7.  The  shoes  or  pads  are  equally  spaced  about  the  shaft 
with  120°  between  pad  centers.  The  three  pad  arrangement  shown  in  Figure  7  is 
mounted  on  a  leveling  washer  which  allows  loads  to  be  distributed  among  the  three 
pads.  The  bearing  inside  diameter  is  3.5  in.  and  the  outside  diameter  is  7  in. 

The  net  area  of  the  three  pads  is  12.3  in. 2  Each  pad  is  fabricated  of  steel,  is 
about  3/4-in.  thick  and  has  a  1/16-in.  thick  babbitt  surface  (ASTM  B23)  and  a 
hardened  steel  insert  in  the  backing  (see  Figure  8).  The  hardened  steel  insert  has 
a  spherical  surface,  which  allows  the  shoe  to  pivot  and  form  a  load-carrying  lubri¬ 
cant  film  between  the  shoe's  babbitt  surface  and  the  runner  surface.  The  two  sets 
of  thrust  pads  were  positioned  so  that  the  pads  from  one  set  were  located  directly 
opposite  the  other  set.  Figure  9  shows  the  bearing  arrangement  with  the  top  load 
ring,  top  set  of  pads,  top  runner  surface,  runner,  bottom  runner  surface,  bottom  set 
of  pads,  and  bottom  load  ring.  Three  load  pistons,  one  over  each  of  the  top  pads, 
apply  a  hydraulic  load  to  the  hearing. 
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8  LEVELING  WASHER 


The  crowning  and  longitudinal  roughness  on  the  bearing  surfaces  were  obtained 
in  the  following  manner: 

•  The  thrust  pads  were  installed  in  their  respective  load  ring  as  though 
ready  to  run. 

•  The  assembly  was  placed  with  the  bearing  surface  down  on  a  flat 
circular  disc  covered  with  emery  paper. 

•  The  disc  was  rotated  while  maintaining  light  pressure  on  assembly. 

The  assembly  was  held  stationary. 

•  Water  was  used  on  the  emery  paper. 

•  The  disc  was  rotated  in  both  directions  until  the  entire  surface  con¬ 
tacted  the  paper,  yielding  a  bright  babbitt  surface. 

•  F.ach  thrust  pad  was  then  measured  to  determine  the  curvatures  and 
roughness . 

TOPOGRAPHICAL  ANALYSIS  OF  THE  TEST  BEARINGS 

The  surface  profiles  and  roughness  of  all  bearing  babbitt  surfaces  and  thrust 
runners  were  obtained  using  the  Talysurf  10  machine  manufactured  by  Rank  Taylor 
Hobson.  It  was  used  without  a  skid  device  so  that  profile  data  could  be  obtained. 
Figure  10  shows  the  standardized  grid  network  set  up  for  the  bearing  and  thrust 
plates.  Profile  traces  were  recorded  for  each  grid  line.  Five  traces  were  taken 
on  each  bearing  pad,  and  four  traces  were  taken  on  each  thrust  runner.  The  trace 
lengths  for  the  thrust  pad  were  approximately  32  and  50  mm  for  the  I  and  0  traces 
and  30  mm  for  the  L,  C,  and  R  traces  (See  Figure  10).  All  runner  traces  were  about 
50  mm.  The  profiles  of  the  surfaces  were  recorded  using  magnification  of  10,000X 
in  the  vertical  direction  and  magnification  of  2,  10,  and  100X  in  the  horizontal 
direction.  Sample  traces  are  shown  in  Figures  11  and  12  for  the  thrust  pad  and 
runner,  respectively.  The  date  and  time,  the  horizontal  and  vertical  magnification 
part  identification,  measurement  location,  and  the  start  and  stop  location  were 
recorded  on  each  trace. 

Then,  the  R^  values  over  the  lengths  of  each  grid  were  recorded.  The  separate 
values  were  recorded  for  each  line  on  the  thrust  runner  and  five  each  on  the  thrust 
pad.  Different  sections  of  the  surface  were  sampled  for  the  different  values. 


Figure  11  -  Thrust  Pad  Talysurf  Traces  at  10,000*  for  Vertical  and  2*  Horizontal 
Magnification  for  Bearing  T-2,  Outside  Axial  Trace 


Figure  12  -  Thrust  Runner  Talysurf  Traces  at  10,000*  Vertical  and  2x  Horizontal 
Magnification  for  Thrust  Runner  A  at  270°  Position 
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Each  Talysurf  10  trace  was  characterized  by  the  following  parameters:  crown 
height,  maximum  peak  to  valley  height,  and  estimated  peak  to  valley  asperity  repre¬ 
sentative  of  surface.  Figure  13  shows  the  basis  on  which  these  parameters  were 
evaluated.  These  values  were  tabulated  for  each  trace  in  addition  to  the  average 
reading,  the  minimum  R^,  and  maximum  R^  readings  for  each.  The  data  for  the  bearing 
pad  (Tables  1  and  2)  and  for  the  thrust  collars  (Table  3)  were  used  to  generate 
input  to  the  computer  program.  The  crown  height  in  the  direction  of  motion  was 
obtained  for  each  thrust  pad  by  taking  the  maximum  value  recorded  for  either  trace 
I  or  0  in  Figure  10.  The  surfaces  of  the  thrust  runner  were  assumed  to  be  flat. 

The  radial  traces  on  the  thrust  shoes  and  runners  were  used  to  determine  the 
asperity  height  magnitude  as  follows:  The  maximum  2c  values  for  the  traces  on  the 
two  runner  surfaces  and  the  thrust  pads  were  obtained  from  Tables  1,  2,  and  3.  The 
resulting  2c  value  for  each  is  shown  in  Table  4.  The  set  of  three  thrust  pads 
identified  with  letter  "T”  were  mated  with  thrust  collar  A;  the  pads  marked  B  were 
paired  with  thrust  collar  B.  A  composite  asperity  height  was  determined  by  averag¬ 
ing  the  particular  pad  and  the  mating  runner.  The  resulting  asperity  height,  2c, 
and  the  crown  height  used  in  the  computer  program  are  shown  in  Table  5. 


MAX  PEAK  TO  VALLEY 


ENTIRE  SURFACE) 


Figure  13  -  Characterization  of  Bearing  Surface  Traces  on  Talysurf  10 
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TABLE  1  -  ROUGHNESS  AND  CROWN  MEASUREMENTS  FOR  BOTTOM 
SET  OF  THRUST  PADS  ON  A  TILT  PAD  BEARING 


Bearing 

No. 

♦  * 

Location 

** 

ra 

RAMin 

R^Max 

Max  *** 

P  to  P 

2c 

Axial 

Crown' 

Radial 
Crown t 

B-l 

R 

12.4 

10.1 

17.9 

430 

30 

- 

-100 

B-l 

C 

13.0 

11.0 

16.8 

150 

40 

- 

-160 

B-l 

L 

21.8 

20.5 

22.9 

190 

50 

- 

-160 

B-l 

I 

- 

- 

- 

240 

45 

+60 

- 

B-l 

0 

- 

- 

- 

165 

50 

+50 

- 

B-2 

R 

13.7 

10.5 

20.0 

130 

30 

- 

-80 

B-2 

C 

11.3 

9.4 

12.5 

130 

35 

- 

-100 

B-2 

L 

15.7 

13.8 

17.9 

190 

50 

- 

-150 

B-2 

I 

- 

- 

- 

200 

45 

+75 

- 

B-2 

0 

- 

- 

- 

350 

40 

-50 

- 

B-3 

R 

16.9 

12.5 

24.0 

200 

50 

- 

-130 

B-3 

C 

15.2 

13.9 

17.5 

150 

40 

- 

-130 

B-3 

L 

30.5 

28.5 

33.0 

600 

45 

- 

+70 

B-3 

I 

- 

- 

- 

450 

40 

+280 

- 

B-3 

0 

- 

- 

- 

310 

60 

+550 

- 

*See  figure  10  for  location  positions 
**Average  of  five  readings  across  surface 


Max  Peak  to  Peak  (P  to  P) 


tPlus  indicates  a  convex  profile;  minus  indicates  a  concave  profile 


‘ring 

Jo. 

Location* 

** 

ra 

RAMin 

RAMax 

Max  *** 

P  to  P 

_  *** 

2c 

Axial 

Crown' 

Radial 

Crownt 

r-i 

R 

10.9 

8.2 

13.3 

200 

20 

- 

+150 

’-i 

C 

16.7 

9.0 

20.0 

550 

40 

- 

+100 

-i 

L 

30.4 

18.0 

73.0 

1150 

60 

- 

+100 

-i 

I 

- 

- 

- 

700 

40 

+160 

- 

’-i 

0 

- 

- 

- 

250 

20 

+200 

- 

-2 

R 

10.7 

9.2 

14.0 

180 

25 

- 

+70 

-2 

C 

11.6 

7.0 

14.5 

300 

25 

- 

+200 

’-2 

L 

16.5 

12.2 

23.5 

650 

40 

- 

+320 

-2 

I 

- 

- 

- 

125 

30 

+80 

- 

-2 

0 

- 

- 

- 

180 

25 

+160 

- 

r-3 

R 

18.1 

13.9 

21.0 

150 

40 

- 

+80 

-3 

C 

18.0 

12.8 

22.0 

180 

35 

- 

+40 

-3 

L 

18.8 

17.3 

20.0 

130 

45 

- 

+20 

-3 

I 

- 

- 

- 

170 

55 

+80 

- 

-3 

0 

- 

- 

- 

170 

25 

+80 

- 
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s  across 

tions 

surface 

0K5S5SS5 


TABLE  3  -  ROUGHNESS  AND  CROWN  OF  THRUST  COLLARS 


TABLE  4  -  MAXIMUM  RECORDED  2c 
VALUES  FOR  THE  THRUST  RUNNER 
AND  BEARING  SURFACES 


Thrust  Runner 
A-l 
B-l 


Maximum  2c 
10 
25 


Bearing  Number 
T-l 
T-2 
T-3 
B-l 
B-2 
B-3 


Maximum  2c 
60 
40 
45 
50 
50 
50 


TABLE  5  -  COMPUTER  INPUT  BASED  ON 
TOPOGRAPHY  OF  SURFACES 


Bearing 

No. 

Crovm 

2c 

T-l 

200 

35 

T-2 

160 

25 

T-3 

80 

27.5 

B-l 

60 

37.5 

B-2 

75 

37.5 

B-3 

550 

37.5 

WWrXWWWWfXn 


COMPUTER  ANALYSIS 


The  program  described  in  the  section  on  theory  was  run  using  the  data  from  the 
topographical  analysis  presented  in  Table  5.  The  following  input  data  concerning 
the  lubricant,  materials,  and  test  conditions  were  used: 

Rotational  shaft  speed  =  35  rpm 

Bearing  mean  diameter  =  (Bore  +  OP)  =5.2  in. 

2 

Lubricant  Viscosity  =  6xl0-^  Reyns 

Shoe  Length  =  2.31  in. 

Asperity  Height  =  value  from  Table  5 

Yield  Stress  =  9000  psi 

Shear/Yield  =  0.2 

Crown  Height  =  value  from  Table  5 

Pivot  Center  =  0.5,  0.32 

Typical  output  for  operation  in  the  hydrodynamic  region  is  shown  in  Figure  14 
and  for  mixed  lubrication  conditions  in  Figure  15.  Two  cases  are  presented  in  each 
figure;  one  allows  negative  pressure  and  the  other  does  not  (imposing  Reynolds- 
Swif t-Stieber  boundar>  conditions).  The  pressure  profiles  for  each  case  are 
presented  in  Figures  16,  17,  18,  and  19. 

The  total  friction  of  the  six  thrust  pads  was  summed  and  then  divided  by  total 
load  to  obtain  the  friction  coefficient  for  a  series  of  VU/PB  values  and  plotted  for 
comparison  with  experimental  results.  Figure  20  shows  the  resulting  friction  curve 
and  Figure  21  the  minimum  film  thickness  predicted  by  the  computer  program  as  a 
function  of  UU/PB  when  negative  pressures  are  allowed  and  not  allowed. 


THIS  CALCULATION  ALLOWED  NEGATIVE  PRESSURES 

SLOPE=l.  74578497  IE-4  PIV0T=0.  5,  0.32 

CENT.  PRESS=0. 499870509989 

PURE  HYDRO  LOAD=1152 .87613215 

MIXED  HYDRO  LOAD=0  PSI 

CONTACT  LOAD=0  PSI 

TOT.  LOAD=1152 .  87613215  PSI 

PURE  HYDRO  FRICT=0 .455753962741  PSI 

MIXED  HYDRO  FRICT=0  PSI 

CONTACT  FRICT=0  PSI 

TOT.  FRICT=0. 455753962741  PSI 

U-ETA/W=2. 146978267E-8 

COEF.  OF  FRICTION  =5 .698976039E-4 

SHOE  SPEED-9.52952.  IN/S 

PRESS.  MOM.  =-0.0663837265664  FT  LB/ IN 

SHEAR  MOM.  =0.0648519658822  FT-LB/lN 

H1/H2  RATI0=3. 70320940581 

H2=l.  491842724 E-4 

C/H2=0.  234609181239 

ABRAMOV IT Z' S  PHI=1.  34062389279 

MIN.  MACRO  FILM  HEIGHT=1 .  0E-41N.  AT  X/B=0 .752047705209 
NO.  ITERATI0NS=4 
NO  MIXED  LUBRICATION 
PLOT  //2 

THIS  CALCULATION  DID  NOT  ALLOW  NEGATIVE  PRESSURES 

SL0PE=1 .564805208E-4  PIVOT  =0.  5,  0.32 

CENT.  PRESS=0. 499974779391 

PURE  HYDRO  LOAD=1206 .30136639 

MIXED  HYDRO  LOAD=0  PSI 

CONTACT  L0AD=0  PSI 

TOT.  LOAD=1206. 30136639  PSI 

PURE  HYDRO  FRICT=0.  450305088924  PSI 

MIXED  HYDRO  FRICT=0  PSI 

CONTACT  FRICT=0  PSI 

TOT.  FRICT=0. 450305088924  PSI 

U-ETA/W-2.  051891898E-8 

COEF.  OF  FRICTION  =5.297745429E-4 

SHOE  SPEED-9.52952  IN/S 

PRESS.  MOM.  —0.0135286387374  FT  LB/ IN 

SHEAR  MOM.  =0.0640766129335  FT-LB/lN 

H1/H2  RATIO-3.  25782683417 

H2-1.  600964244E-4 

C/H2-0. 218618249131 

ABRAMOV ITZ' S  PHI-1.  24924713789 

MIN.  MACRO  FILM  HEIGHT-1.  OE-4IN.  AT  X/B-0.  725918751951 
NO.  ITERATIONS- 3 

REYNOLDS-SWIFT-STIEBER  CAVITATION  CONDITION  IMPOSED 
NO  MIXED  LUBRICATION 
PLOT  #3 


Figure  14  -  Typical  Computer  Output  Under 
Hydr idynamic  Lubrication  Conditions 


THIS  CALCULATION  DID  NOT  ALLOW  NEGATIVE  PRESSURES 

SLOPE=l.  15566044E-4  PIVOT  =0.  5,  0.32 

CENT.  PRESS=0 •  499920907264 

PURE  HYDRO  LOAD=3679.  03496371 

MIXED  HYDRO  LOAD=1386.  29530187  PSI 

CONTACT  LOAD=0. 19443687 2283  PSI 

TOT.  L0AD=5065.  52470245  PSI 

PURE  HYDRO  FRICT=0.  9356570736  PSI 

MIXED  HYDRO  FRICT=0.  246750572852  PSI 

CONTACT  FRICT=0 .  0388873744566  PSI 

TOT.  FRICT=1 .  22129502091  PSI 

U-ETA/W=4.  886364484E-9 

COEF.  OF  FRICTI0N=3 .566654547E-4 

SHOE  SPEED=9 .  52952  IN/S 

PRESS.  MOM.  =-0.178157352734  FT  LB/IN 

SHEAR  MOM.  =0.173785396295  FT-LB/lN 

H1/H2  RATI0=3 .  24726979178 

H2=l .  187919504E-4 

C/H2=0.  294632758297 

ABRAMOV ITZ' S  PHI=1.  6836157617 

MIN.  MACRO  FILM  HEIGHT=3 .  0E-5IN.  AT  X/B=0.  666848475978 
NO.  ITERATI0NS=4 

REYNOLDS-SWIFT-STIEBER  CAVITATION  CONDITION  IMPOSED 
MIXED  LUBRICATION  WITHIN59AND74 .501% 

plot  no 

THIS  CALCULATION  ALLOWED  NEGATIVE  PRESSURES 

SLOPE=l .  869157 56E-4  PIVOT  =0.5,  0.32 

CENT.  PRESS=0 .  499818277712 

PURE  HYDRO  LOAD=3127.  8232727 

MIXED  HYDRO  LOAD=1254.  44868354  PSI 

CONTACT  LOAD=3 .  72204560221  PSI 

TOT.  LOAD=4385 .  99400184  PSI 

PURE  HYDRO  FRICT=1.  06370550783  PSI 

MIXED  HYDRO  FRICT=0.  465346373011  PSI 

CONTACT  FRICT=0. 744409120441  PSI 

TOT.  FRICT=2 .  27346100128  PSI 

U-ETA/W=5.  64341857-9 

COEF.  OF  FRICTION  =7 .052614264E-4 

SHOE  SPEED=9.  52952  IN/S 

PRESS.  MOM.  =0.354420588933  FT  LB/IN 

SHEAR  MOM.  =0.  323504406638  FT-LB/lN 

H1/H2  RATIO=7.  40885652335 

H2=6 .  737167461E-5 

C/H2=0.  519506160446 

ABRAMOVITZ'S  PHI=2.  96860663112 

MIN.  MACRO  FILM  HEIGHT=2.  5E-5IN.  AT  X/B=0.  7698596227 
NO.  I TE RATION S=4 

MIXED  LUBRICATION  WITHIN66 AND88 .  001% 

PLOT  #11 


Figure  15  -  Typical  Computer  Output  Under 
Mixed  Lubrication  Conditions 


PRESSURE  PSI 


PRESSURE  VS  PAD  WIDTH 


Figure  16  -  Pressure  Distribution  under  Hydrodynamic 
Lubrication  Conditions  allowing  Negative  Pressures 


PRESSURE  VS  PAD  WIDTH 


Figure  17  -  Pressure  Distribution  under  Hydrodynamic 
Lubrication  Conditions  with  Reynolds-Swif t-Stieber 
Boundary  Conditions 


EXPERIMENTAL  RESULTS 

After  topographical  analysis  of  the  bearing  and  mating  surfaces,  the  bearing 
test  machine  was  assembled.  A  standard  Navy  lubricant  for  thrust  bearings,  2190 
TEP,  was  used  in  the  first  test  sequence.  The  viscosity  at  two  reference  tempera¬ 
tures  (104°F  and  212°F)  and  the  specific  gravity  at  60°F  were  determined  for  a 
sample  from  the  oil  used  in  the  experiment.  The  viscosity  at  other  temperatures 
was  obtained  using  ASTM  standard  viscosity-temperature  charts®  for  liquid  petroleum 
(see  Figure  22).  The  first  runs  were  made  at  ambient  oil  temperatures.  (No 
external  heat  was  added  to  the  lubricant  during  the  test).  The  speed  of  the  test 
machine  was  set  to  a  nominal  35  rpm  with  no  applied  bearing  load.  Loads  were 
applied  in  increments  of  approximately  75  psi  bearing  unit  load  up  to  about  450 
psi.  Torque  and  other  data  were  recorded  immediately  after  a  change  in  load  and 
after  15  min  operation  at  each  load.  Other  parameters  recorded  were  time  of  test, 
rotational  speed,  load  cell  readings,  temperature  of  the  bulk  oil,  temperature  of 
oil  near  bearing  outlet,  water  sump  temperature,  and  bearing  torque.  Table  6  shows 
results  of  the  ambient  temperature  runs. 

The  model  predicts  that  there  should  be  hydrodynamic  behavior  when  PU/PB 
exceeds  2x10"®,  according  to  Figures  20  and  21.  The  lowest  pu/PB  value  reached 
in  the  ambient  test  sequence  was  2xl0-^.  The  lubricant  was  heated  to  150°F  to 
reduce  the  viscosity  so  that  lower  values  of  PU/PB  could  be  reached  while  speed 
and  load  conditions  were  maintained  within  reasonable  limits.  The  lubricant  used 
as  in  the  ambient  temperature  runs  was  heated  by  circulating  water  through  the 
coils  in  the  oil  sump.  Table  7  shows  results  of  the  higher  temperature  runs. 
Although  PU/PB  of  3.5x10"®  was  reached,  the  friction  data  still  appeared  to  be 
obeying  hydrodynamic  principles. 


Test 

Descri pt i on 

Time 

On 

(hr) 

Time 

Off 

(hr) 

Speed 

N  U 

(rpm)  in/s 

Load  Cell 
Readings 
(lb) 

Total 

Load 

Cells 

(lb) 

Bearing 
Load 
<  psl) 

Bulk  Oil 
Tempera¬ 
ture 

Cf) 

Oil 
Vis¬ 
cosity 
( cs) 

Oil 

Viscosity 

(Reyns) 

UU/PB 

Torque 

in-lb 

* 

f=— L- 
2  FR 

3  Pad  Thrust 

1310 

35 

9.63 

333 

280 

320 

933 

75.85 

73 

260 

32X10-6 

1.74X10'6 

72.23 

14.7X10*3 

Bearing 

1323 

36 

9 . 5" 

in  3 

250 

290 

843 

68.54 

74 

260 

32X10-& 

2X10-6 

63.89 

14.4X10')  | 

with  2190 

TEP  Oil 

1328 

38 

9.0. 

►•10 

520 

590 

1720 

139.84 

75 

250 

31X10-6 

9 .0X10-7 

108.34 

12X10-? 

at  Room 

1343 

36 

9.9  • 

585 

500 

560 

1645 

133.74 

76 

250 

31X10-6 

1.0X10-6 

97.23 

11.3X10"3 

Tempera- 

ture 

1347 

36 

9.90 

950 

820 

910 

2680 

217.9 

76 

240 

30X10-6 

5.8X10"7 

136.12 

9.7X10-5 

1402 

36 

9.90 

930 

800 

880 

2610 

.12.2 

77 

240 

30X10-6 

5.99X10-7 

125.01 

9.1X10-3 

1403 

36 

9.90 

1250 

1090 

1200 

3540 

287.8 

77 

220 

27X10-6 

3.98X10*7 

144.46 

7 .8X10-5 

14  20 

36 

9.90 

1230 

1080 

1180 

3490 

283.7 

79 

220 

27X10-6 

4.03X10*7 

130.57 

7.1X10-5 

1423 

36 

9.90 

1570 

1380 

1520 

4470 

363.4 

79 

210 

26X10-6 

3.03X10-7 

161.12 

6.8X10-5 

1435 

36 

9.90 

1560 

1375 

1500 

44  35 

360.6 

80 

210 

26X10-6 

3.06X10'7 

144.46 

6.2X10-5 

1440 

36 

9.90 

2040 

1810 

1990 

5840 

474  .8 

81 

200 

24X10-6 

2.14X10-7 

183.35 

6.0X10-5 

1450 

36 

9.90 

2040 

1800 

1980 

5820 

473.2 

82 

24X10-6 

2.15X10-7 

161.12 

5.3X10-5 

TABLE  7  -  DATA  FROM  RUN  2  -  2190  TEP  AT  150°F 


Test 

Descript  ion 

Time 

On 

(hr) 

Time 

Off 

(hr) 

Speed 

N  t! 

(rpm)  in/s 

Load  Cell 
Readings 
(lb) 

Total 

Load 

Cells 

(lb) 

Bulk  Oil 
Bearing  Tempera- 
Load  ture 

(PSD _ CF) 

Oil 

Vis¬ 

cosity 

(cs) 

Oil 

Viscosity 

(Revns) 

UU/PB 

Torque 

in-lb 

2FR 

3  Pad  Thrust 

1257 

36.5 

10.04 

340 

300 

292 

932 

75.77 

s^xio-p 

1.87X10"' 

16.67 

3.41X10-7 

Bearing 

1312 

36.5 

10.04 

308 

267 

255 

830 

67.48 

153 

25 

3.1X10-6 

1.97X10"7 

12.22 

2.81X10-3 

with  2190 

TEP  Oil 

1316 

36 

9.90 

610 

550 

545 

1705 

138.62 

154 

24 

2 .9X10*6 

8.87X10'8 

27.78 

3.10X10-3 

at  1 50°F 

1331 

36 

9.90 

590 

520 

520 

1630 

132.52 

155 

24 

2.9X10-6 

9 .27X10*6 

22.22 

2.6X10*3 

1335 

36 

9.90 

935 

840 

845 

2620 

213 

154 

24 

2 .9X10*6 

5.77X10*8 

38.89 

2.83X10-3 

1350 

36 

9.90 

890 

800 

800 

2490 

202.4 

145 

30 

3.6X10’6 

7.54X10*6 

36.11 

2.76X10*3 

1354 

36 

9.90 

1240 

1125 

1150 

3515 

285.8 

144 

30 

3.6X10-6 

5.34X10'8 

50 

2.71X10*3 

1409 

36 

9.90 

1190 

1080 

1110 

3380 

274.8 

143 

30 

3.6X10-6 

5.55X10'8 

44.45 

2.50X10-3 

1412 

36 

9.90 

1545 

1400 

1445 

4  390 

357 

143 

30 

3.6X10-6 

4.27X10"8 

58.34 

2.53X10*3 

1427 

36 

9.90 

1525 

1380 

1440 

4  345 

353.3 

146 

30 

3.6X10*6 

4.32X10*6 

50 

2.19X10-3 

1429 

9.90 

1875 

1690 

1785 

5350 

4  35 

147 

28 

3.5X10-6 

3.41X10*6 

63.89 

2.27X10*3 

1444 

36 

9.90 

1870 

1710 

1775 

5355 

435.4 

150 

26 

3.2X10-6 

3.11X10"8 

50 

1.78X10-3 

34 


A  low  viscosity  fluid  was  sought  to  extend  the  test  to  lower  values  of  PU/PB. 

A  hydraulic  fluid  meeting  MIL-H-5606^  was  selected.  The  viscosity-temperature 
relationship  for  a  sample  of  the  MIL-H-5606  fluid  was  obtained  and  plotted  in 
Figure  22.  These  fluids  exhibit  non-Newtonian  behavior  and  break,  down  under  high 
shear  rates. *0  While  the  fluid  is  non-Newtonian  at  about  59°F,  its  behavior  at 
100°F  approaches  Newtonian.  The  shear  rate  at  pU/PB  of  2xl0~^  for  the  computerized 
model  with  Reynolds ' -Swif t-Stieber  cavitation  condition  p:oduc-js  a  shear  rate  of 
16xl04  reciprocal  seconds.  The  viscosity  loss  due  to  shear  should  amount  to  only 
10%.  The  fluid  was  heated  to  about  100°F  and  tests  conducted  in  the  bearing  machine 
ir  a  similar  manner  to  those  described  earlier.  The  results  of  those  tests  are 
listed  in  Table  8. 

The  data  from  the  three  test  runs  are  compared  to  the  computer-generated  results 
in  Figure  25.  The  friction  coefficient  is  defined  by  the  following  expression: 

T  =  f  FR 

where:  T  is  the  torque  in  in-lb, 

f  is  the  friction  coefficient, 

F  is  the  total  applied  thrust  load, 
and  R  is  the  radius  to  the  pad  center. 

The  measured  torque  on  the  test  machine  is  divided  by  two  because  it  is  the  result 
of  two  sets  of  bearings.  The  data  for  the  2190  TEP  oil  in  both  ambient  and  150°F 
runs  are  parallel  to  the  predicted  values  and  apparently  operating  hydrodynamically , 
but  the  magnitude  of  the  friction  measured  is  about  four  times  the  value  predicted. 
The  results  with  the  MIL-H-5606  fluid  indicate  that  the  bearings  were  operating  in 
mixed  lubrication  conditions. 


TABLE  8  -  DATA  FROM  RUN  3  -  MIL-H-5606  FLUID  AT  105  F 


Test 

Descript  ion 

Time 

On 

W 

T  ime 
Off 
(hr) 

Speed 

N  U 

(rpm)  in/s 

Load  Cell 
Readings 
(lb) 

Total 

Load 

Cells 

(lb) 

Bearing 

Load 

(psl) 

Bulk  Oil 
Tempera¬ 
ture 

(  °F) 

Oil 

Vis¬ 

cosity 

(cs) 

Oil 

Viscosity 
( Reyns) 

uU/PB 

Torque 

in-lb 

* 

f«-I_ 

2FR 

3  Pad  Thrust 

914 

32 

8 .8 

370 

335  325 

1030 

83.75 

14 

1.7X10-6 

7.64X10-8 

37.78 

7X10"-* 

Bearing 

929 

30 

8.25 

320 

285  270 

875 

71.14 

105 

14 

1.7X10-6 

8.44X10-8 

15 

3.3X10"3 

with  MIL- 

H-5606 

934 

29 

7.98 

625 

570  560 

1755 

142.7 

106 

14 

1.7X10-6 

4.07X10"8 

96.12 

10.4X10"3 

Fluid 

949 

29.5 

8.11 

585 

530  525 

1640 

133.3 

107 

13 

1.6X10"6 

4 .17X10”° 

62.78 

7.3X10_3 

933 

29.5 

8.11 

930 

850  870 

2650 

215.5 

107 

13 

1.6X10-6 

2.10X10-8 

165.6 

11.9X10-3 

36 


DISCUSSION 


The  friction  coefficients  measured  under  hydrodynamic  conditions  during  the 
experiment  using  the  2190  TEP  turbine  oil  were  about  four  times  those  predicted  by 
the  computer  model  under  hydrodynamic  conditions.  Both  the  ambient  and  150°F 
temperature  runs  using  that  oil  exhibited  friction  representative  of  a  complete 
hydrodynamic  lubricant  film  separating  the  surfaces. 

We  suspect  that  of  che  assumptions  made  in  the  theoretical  model,  that  for 
neglecting  of  side  leakage  would  have  influenced  the  results  most  significantly. 
Fuller^!  pointed  out  the  significant  influence  of  side  leakage.  He  presented  side 
leakage  correction  factors  for  the  load  and  film  thickness  expressions  describing 
a  simple  fixed  inclined  flat  bearing  that  has  an  inlet  film  thickness  of  double 
the  outlet  film  thickness.  He  modified  his  derived  expressions  for  bearing  load 
and  minimum  film  thickness  by  a  side  leakage  factor,  n,  as  follows: 


6uUl  cPK. 


and  h  = 
o 


bwuinKp 

p  _ 


The  side  leakage  factors  were  based  upon  experiments  of  Kingsbury  and  Needs 
(as  cited  by  Fuller^)  and  are  a  function  of  the  bearing  length  to  width  ratio. 
The  coefficient  of  friction  for  Fuller's  bearing  configuration  is  given  by: 


yBlUKf 


f  =  —  = 

W 


It  is  proportional  to  n According  to  Fuller,  the  correction  factor,  n,  for 
side  leakage  for  the  thrust  bearings  evaluated  in  our  work  is  0.32.  Thus  the  cal¬ 
culated  friction  should  be  multiplied  by  (0.32)”^/2  or  1,77  to  yield  a  better 
value  for  comparison  with  measured  friction. 

Correcting  for  side  leakage  would  account  for  nearly  half  of  the  difference 
between  our  predicted  and  measured  results.  In  addition,  the  thrust  pads  used  in 
our  study  were  crowned  both  in  the  direction  of  motion  and  in  the  radial  direction. 
Under  these  conditions,  it  is  conceivable  that  side  leakage  correction  could 


account  for  the  four-fold  difference  between  the  friction  coefficients. 


The  model  assumed  a  parabolic  surface  for  each  bearing.  The  topographical 
measurements,  however,  showed  that  not  all  bearing  surfaces  could  be  represented 
by  a  simple  parabola.  The  capacity  to  input  varying  shapes  of  the  bearing  surface 
is  needed.  Accordingly  we  have  begun  work  with  Dr.  C.  Pan  of  Columbia  University 
(N00167-82-C-0134)  which  includes  studies  of  side  leakage  and  variable  pad  shapes. 
Dr.  Pan's  model  is  a  thrust  pad  with  a  20  x  20  finite  element  grid.  Surface 
shapes  can  be  set  by  the  initial  values  given  to  the  400  grid  points. 

Thermo-elastic  effects  should  be  considered  in  the  development  of  a  more 
sophisticated  model  because  thermal  and  elastic  distortions  change  the  shapes  of 
the  thrust  pad  surfaces.  Some  work  has  been  done  in  these  areas.  A  finite  element 
model  has  been  developed  to  determine  the  distortions  due  to  thermal  gradients 
within  the  thrust  pad. 12  under  normal  (hydrodynamic)  conditions,  lubricant  shear 
and  minor  asperity  contact  generate  heat.  Some  of  the  heat  is  absorbed  by  the 
fluid  or  through  the  load  zone;  the  rest  is  carried  away  through  the  thrust  pad. 

The  surface  temperature  at  the  film  would  be  hottest  and  there  is  a  thermal  gra¬ 
dient  through  the  thickness  of  the  pad.  This  gradient  would  Increase  the  initial 
crown  and  the  crown  caused  by  pressure  loading  on  the  face.  Thus,  heat  transfer 
must  be  analyzed  to  develop  the  temperature  profile  in  the  thrust  pad.  A  pre¬ 
liminary  finite  element  study  of  heat  transfer  under  hydrodynamic  lubrication^  has 
been  conducted  under  this  program.  A  general  analysis  of  interactive  mechanisms 
of  sliding-surface  bearings  was  presented  by  Daugherty  and  Pan.^  However, 
additional  work  is  needed  in  elasto-hydrodynamics ,  before  the  results  can  be 
incorporated  into  our  model. 

The  tests  conducted  with  the  MIL-H-5606  fluid  resulted  in  frictional  behavior 
typical  of  the  mixed-lubrication  regime.  The  non-Newtonian  effects  were  much 
greater  than  expected.  In  Figure  23,  the  friction  coefficients  were  plotted  against 
a  wU/PB  value  based  upon  measurements  of  the  kinematic  viscosity  of  the  fluid.  The 
kinematic  viscosity  is  determined  by  measuring  the  time  for  the  liquid  to  flow  under 
gravity  through  a  glass  viscometer  and  thus  does  not  consider  high  shear  rates  such 
as  those  of  the  experiment.  The  viscosity  of  the  MIL-H-5606  fluid  would  need  to  be 
about  one  third  of  the  value  used  for  the  friction  data  shown  in  Figure  23  to  fit 
the  curve  expected. 

Asperity  interactions  occurred  during  the  runs  with  MIL-H-5606  fluid.  Wear 
was  apparent  during  the  15  min  operation  at  a  given  load.  The  resulting  decrease 
in  friction  coefficient  may  have  been  caused  by  the  polishing  of  high  asperities. 


The  arrows  in  Figure  23  show  the  sequence  of  tests  conducted  for  the  MIL-H-5606 
fluid.  The  friction  coefficient  decreased  after  15  min  under  given  load,  then 
increased  when  additional  load  was  applied.  Furthermore,  the  model  as  presently 
structured  does  not  consider  the  contribution  of  the  mechanical-chemical  inter¬ 
actions  of  the  surfaces  with  the  fluid  and  its  additives  that  occur  in  the  mixed 
and  boundary  lubrication  zones. 

Studies  are  being  conducted  at  Mississippi  State  University  by  Dr.  A.  Stiffler 
(DTNSRDC  Contract  N61533-83-M-2799)  to  predict  the  load  support  and  friction  coef¬ 
ficient  in  sliding  pairs  when  an  ensemble  of  asperities  is  operating  in  near 
contact.  Dr.  Stiffler's  analysis  will  include  one  dimensional  microtopographic 
studies  for  an  asperity  against  a  flat  surface  and  an  asperity  against  another 
asperity,  for  elastic  deformations  of  asperities,  for  non-Newtonian  lubricants, 
and  for  temperature  effects. 

CONCLUSIONS 

Measured  friction  of  a  crowned  tilt-pad  thrust  bearing  with  longitudinal 
roughness  when  plotted  against  the  uU/PB  parameter  was  parallel  to  but  higher  than 
the  results  predicted  by  the  model  under  hydrodynamic  lubrication  conditions. 

We  now  believe  that  major  differences  were  due  primarily  to  the  neglect  of  side 
leakage  in  the  model.  Additional  improvements  could  be  made  to  the  model  by  consid¬ 
ering  variable  viscosity  in  the  film  and  different  film  shapes.  The  present  model 
cannot  handle  the  complex  tribological  interactions  that  occur  in  experiments  in  the 
mixed  and  probably  in  the  boundary  lubrication  regimes  of  tilt-pad  thrust  bearings. 

RECOMMENDATIONS 

The  present  model  should  be  expanded  to  consider  side  leakage  effects,  variable 
viscosity,  thermo-elastic  deformations,  pad  shapes,  skewed  roughness  patterns,  and 
improved  local  asperity  configurations  for  a  tilt-pad  bearing  operating  under  hydro- 
dynamic  and  mixed  lubrication.  Also,  more  detailed  studies  are  needed  of  the  local 
behavior  of  asperities  under  concentrated  loads.  Finally,  additional  experiments 
are  needed  to  verify  some  of  the  assumptions  in  the  present  model.  Such  experiments 
will  isolate  the  effects  of  selected  parameters  on  the  tribological  system  and  lead 
to  a  better  understanding  of  the  mechanisms  involved,  be  they  physical,  chemical  or 
(as  expected)  a  combination  of  both. 
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